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ABSTRACT 

Context. Simulations of astrophysical disks in the shearing box that are subject to the magnetorotational instability 
(MRI) show that activity appears to be reduced as the magnetic Prandtl number P m is lowered. It is therefore important 
to understand the reasons for this trend, especially if this trend is shown to continue when higher resolution calculations 
are performed in the near future. Calculations for laboratory experiments show that saturation is achieved through 
modification of the background shear for P m <C 1. 

Aims. Guided by the results of calculations appropriate for laboratory experiments when P m is very low, the stability of 
inviscid disturbances in a shearing box model immersed in a constant vertical background magnetic field is considered 
£\j j under a variety of shear profiles and boundary conditions in order to evaluate the hypothesis that modifications of the 

shear bring about saturation of the instability. Shear profiles q are given by the local background Keplerian mean, qo, 
plus time- independent departures, Q(x), with zero average on a given scale. 

Methods. The axisymmetric linear stability of inviscid magnetohydrodynamic normal modes in the shearing box is 
j analyzed. 

Results, (i) The stability/instability of modes subject to modified shear profiles may be interpreted by a generalized 
i-C ■ Velikhov criterion given by an effective shear and radial wavenumber that are defined by the radial structure of the 

Oh' mode and the form of Q. (ii) Where channel modes occur, comparisons against marginally unstable disturbance in the 

classical case, Q = 0, shows that all modifications of the shear examined here enhance mode instability, (iii) For models 
with boundary conditions mimicing laboratory experiments, modified shear profiles exist that stabilize a marginally 
unstable MRI for Q = 0. (iv) Localized normal modes on domains of infinite radial extent characterized by either single 
defects or symmetric top-hat profiles for Q are also investigated. If the regions of modified shear are less (greater) than 
the local Keplerian background, then there are (are no) normal modes leading to the MRI. 

Conclusions. The emergence and stability of the MRI is sensitive to the boundary conditions adopted. Channel modes do 
not appear to be stabilized through modifications of the background shear whose average remains Keplerian. However, 
systems that have non-penetrative boundaries can saturate the MRI through modification of the background shear. 
\Q ' Conceptually equating the qualitative results from laboratory experiments to the conditions in a disk may therefore be 

< ^* ) misleading. 
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1. Introduction depend on P m : as this quantity decreases the vigor of the 

turbulent state decreases and if P m drops below some criti- 

. . ; After almost 40 years of investigation the source of the cal value> then the turbu i en ce appears to vanish altogether. 

anomalous transport in accretion disks still remains an open 
k> ' question. As it is mostly assumed that the transport is a Estimates of the properties of cold astrophysical disks, 
; reflection of some underlying turbulent state, identification hke protoplanetary disks, show that their characteristic 
. of its source mechanism and process has been the focus of p m' s are 011 tne order of 10 at best and that the MRI 
- - ' much research activity. A leading candidate mechanism is ma y be operating only in the disk's corona as its ionization 
the magnetorotational instability (Balbus & Hawley 1991). fraction is sufficient to merit an MHD description there 
The magnetic Prandtl number (P m ), i.e. the ratio of a rather than near the disk midplane (see Balbus & Henri 
fluid's viscosity to its magnetic diffusivity, appears to play 2008 and references therein). Thus, the fate of the MRI m 
an important role in the ability of the MRI to drive a fully- disk environments m which the magnetic Prandtl number 
developed turbulent state. Recently resolved simulations 18 characteristically small remains to be settled, 
of the MRI in a shearing box environment (e.g. Lesur & It is not clear why turbulent transport appears to van- 
Longaretti 2007; Fromang et al. 2008) show the appearance ish in numerical experiments when P m weakens. One pos- 
of a turbulent state for moderately high Reynolds numbers sibility may be that under those conditions the MRI can- 
(Re) when P m is an order 1 quantity or higher. However, not grow sufficiently to excite a secondary transition that 
the amount of transport delivered by the MRI appears to would, in turn, open the way to a turbulent cascade. 
Theoretical considerations of the MRI for restrictive con- 
Send offprint requests to: O.M. Umurhan e-mail: figurations like laboratory experiments with cylindrical ge- 
umurhan@maths . qmul .ac.uk ometry predict that a low P m sheared fluid undergoing the 
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MRI settles onto a pattern state whose momentum trans- 
port scales as ~ (P m ) a where a > 1 - in other words, 
saturation of the instability is sensitive to the microscopic 
viscosity of the fluid. This counterintuitive result is ratio- 
nalized by noting that a weak fluid viscosity means that it 
takes very little effort to readjust the underlying flow profile 
and establish a new shear in which the instability cannot 
operate. The stronger the viscosity, then, it becomes more 
difficult for the fluid to shut off the shear and the instability 
may operate unabatedly during the fluid's nonlinear devel- 
opment and eventual cascade to turbulence via secondary 
instabilities (e.g. parasitic instabilities, Xu & Goodman, 
1994). 

In the idealized quasi-linear study of Knobloch & Julien 
(2006) it was shown that in the limit of very large hydro- 
dynamic and magnetic Reynolds numbers (i.e. Re 3> 1 
and Re m 3> 1 respectively) the azimuthal velocity pro- 
file tips over to shut off the instability by establishing a 
new velocity profile which reduces the shear throughout 
the domain. Their analysis shows that in the limits where 
Re — > oo and Re m — > oo the resulting net velocity pro- 
file tends toward zero shear. Examinations of the thin-gap 
Taylor-Couette system in the low P m limit (Umurhan et 
al. 2007a; Umurhan et al. 2007b) show that once the sys- 
tem reaches saturation, the resulting velocity field within 
large mid-portions of the experimental domain is character- 
ized by weakened shear. However these regions, where the 
tendency toward instability is reduced, are sandwiched by 
regions of strengthened shear, where the tendency toward 
instability is enhanced. Nonetheless, the aggregate config- 
uration results in a profile that is stable against the MRI 
or any secondary instability and the system settles onto 
non-turbulent pattern state. Equally relevant is that in this 
pattern state the amplitudes of all fluid quantities, except 
the modified shear, scale as some power of the P m . By com- 
parison the modification of the shear is an order 1 quan- 
tity. Thus, in the limit where P m becomes very small only 
the modification to the background shear appears as a no- 
ticeable response in the idealized studies of the laboratory 
setup. 

Because resolved shearing box numerical experiments 
of high Re and low P m flows appropriate for cold astro- 
physical disks is currently out of reach, speculation at this 
stage is justifiable. In particular, what may be the ef- 
fective azimuthal velocity profile of high Re and low P m 
disks (or their sections)? As Ebrahimi et al. (2009) recently 
note, quasi-linear saturation of the shear profile might be 
an unreasonable expectation given that the strong gravity 
tends to restore the shear profile within disks (cf. Julien 
& Knobloch, review in preparation 2009). However, the 
theoretical results of the thin-gap Taylor-Couette system 
and the quasi-linear studies, together with the objection of 
Ebrahimi et al., point to a possible hybrid scenario: What if 
the velocity profile of a low P m disk is on average Keplerian 
but locally exhibits alternating zones of very weak (or no) 
shear and very strong shear? This hypothetical configura- 
tion of the azimuthal velocity profile, an example of which 
is depicted in Figure [TJ might very well be driven into place 
by the MRI under very low P m situations. The recent stud- 
ies of Julien & Knobloch (2006) and Jamroz et al. (2008) 
argue for the relevance of this scenario. 

This study is an examination of the axisymmetric lin- 
ear stability of incompressible ideal magnetohydrodynamic 
disturbances in a shearing box threaded by a constant ver- 




Radial Distance: r 



Fig. 1. Qualitative depiction of a hypothesized mean rota- 
tional velocity profile for a low magnetic Prandtl number 
disk suggested by Jamroz et al. (2008). The dashed line 
shows the usual Keplerian velocity profile describing a ro- 
tationally supported disk. The mean velocity profile result- 
ing from the saturation of the MRI in the low P m limit is 
speculated to consist of alternating regions of low and high 
shear (black lines). The mean shear over the extent of the 
disk follows the general Keplerian profile. 



tical magnetic field and characterized by a variety of mean 
velocity profiles, including an extreme instance of the one 
appearing in Figure [T] Q Such configurations might be rep- 
resentative of real astrophysical disks characterized by very 
small values of P m . 

Guided by the low P m results of Knobloch & Julien 
(2005) and Umurhan et al. (2007a-b), the analysis under- 
taken in this study is inviscid. The only low P m effect that is 
included here is the possibility that the fundamental shear 
profile is significantly modified. This altered shear profile 
may be easily included into the equations of motion as any 
radially dependent barotropic steady shear profile is a per- 
mitted equilibrium flow solution of the small shearing-box 
equations. 

The range of questions that figure prominently are: (i) 
What happens to the unstable mode leading to the MRI 
when there are departures of the steady shear from the 
background Keplerian state? (iii) How does the nature or 
existence of unstable modes depend upon the radial bound- 
ary conditions of the shearing box (e.g. whether it is open 
or periodic in the radial direction)? (ii) Can one infer a 
possible reason or explanation for why turbulence seems to 
vanish in numerical experiments of the shearing box when 
P m is small? 

This work is organized as follows: In Section 2 the equa- 
tions of motion are laid out and the framework within which 
the stability analysis is discussed. Section 3 develops the 
profiles for the steady state for given modifications of the 
background shear while in Section 4 the normal mode sta- 
bility analysis is developed and the governing ODE is de- 
rived including a statement of the generalized Velikhov cri- 
terion (Velikhov 1959). Section 5 presents the results of the 
stability analysis for a variety of boundary conditions and 
shear profiles and their forms. The implications of the re- 
sults are discussed in the final section. 



Extreme in the sense that regions of sharp shear are repre- 
sented by delta-functions in the analysis herein. See also discus- 
sion in Section 6. 
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2. Equations of Motion 

The small shearing box equations (SSB, Lesur & Longaretti 
2007; Regev & Umurhan 2008) are adopted which are the 
incompressible incarnation of the usual shearing box equa- 
tions (Goldreich & Lynden-Bell 1965). Their usage is justi- 
fied as compressibility is an inessential feature for the MRI. 
The box is moving in a rotating frame moving with the lo- 
cal Keplerian velocity V^H The length scales of the box 
are measured by a size L which is much smaller than the 
radial scale of the disk Rq. On these scales, the background 
Keplerian shear appears as a linear Couette flow profile. 
Time units are scaled by the local rotation time of the disk 
CI(Rq), as measured at the fiducial radius Rq,_ implying that 
the velocities of interest are scaled by Uq — LCl(Ro). Given 
this, the disparity of the ratio of the local disk soundspeed 
to the Keplerian speed (since the gas is relatively cold), and 
the assumption that the scale L is much less than the ver- 
tical scale height of the disk, the flow dynamics have the 
character of incompressible rotating magnetic Couette flow. 
Magnetic fields are scaled by a reference field scale B and 
the density is p. These quantities form a velocity defined 
to be the Alfven speed via U 2 
Cowling number is defined 

u 2 A 

c = — 

u 21 

which may be considered as related to the inverse-square 
of the "/3-parameter" of plasma physics. This is not to be 
confused with the parameter /3 used later on in this work 
to designate radial wavenumbers. The non-dimensionalized 
SSB equations of motion are 



B 2 1 Att p. Furthermore the 



(d t + u • V)u + 2f2 z x u = 

— Vp + 
P 

V • u = 

d t B = 

J = 



2nlq xx + C J x B, (1) 

0, (2) 

VxuxB, (3) 

V x B, (4) 



in which u is the velocity vector, p is the pressure, p is the 
constant density (which is 1 in these units), fio is the local 
rotation rate but in the units used here it is also 1. The 
parameter qo defines the amplitude of the local Keplerian 
shear and is formally given by 



R on 
7i ' dR 



Rn 



which is equal to 3/2 for rotationally supported Keplerian 
disks. This is the value assumed for qo throughout the re- 
mainder of this work. The magnetic field is given by the 
vector B and the current J and is governed by Ampere's 
Law.0 

3. Steady State 

It is assumed that the fluid is of constant density po and 
the fluid is threaded uniformly throughout by a constant 



2 Overbars over quantities refer to their dimensional values. 

3 Note that by its construction © contains the statement that 
the time evolution of the divergence of the magnetic field is zero. 
If initially V - B = everywhere, then it remains identically zero 
subsequently. 



vertical magnetic field, i.e. B = B z i. The current source of 
this background field is taken to be external to the domain. 
As stated in the Introduction the steady velocity field that 
represent departures to the Keplerian shear are taken as 
given and represented by U(x). The total mean velocity is 
written as V(x) which is the sum of the Keplerian portion 
—qoflx and the departures U. Thus the steady state con- 
figuration requires satisfying the radial momentum balance 
which becomes 



2f2o(— qo^o x + U) = —d x TiJo + 257q(/o^, 



(5) 



in which the total mechanical and magnetic pressure is 
given by 

= Po/Pa + \CB 2 Z , 

where po is the steady state mechanical pressure field. For 
the model analyzed here where B z is constant, the steady 
pressure field po is related to the azimuthal velocity depar- 
tures by 



d x p = 2Q Up. 



(6) 



The above expression represents a radially geostrophic 
state. Note that the solutions to po are smooth and well 
behaved so long as V is continuous across the domain. In 
the axisymmetric linear stability analysis that follows only 
the steady state total vorticity field, denoted by q, plays a 
role and it is defined as 



Q = Qo + Q, 



(7) 



where Q(x) = —d x U/Clo- A variety of functional forms for 
Q will be examined. The only restriction on Q will be that 
its domain integral is bounded, i.e., 



Qdx 



p 



< oo, 



where T> symbolizes the domain under consideration. In 
many cases this integral will be zero. 

4. Linear Theory 

Axisymmetric infinitesimal disturbances about the steady 
state shaped by q(x) are introduced into the governing 
equations of motion (1-4) revealing, 



dtu — 2VLqv' 
d t v' + (2 - q)n u' 
d t w' 
d x u + d z w 

d t b' x 
d t b' v 
d t b' z 

in which u' i 



~d x (p'/ Pa +CB z b' z ] 
CB z d z b' y , 

-d z (p'/p +CB z b' z ] 
0, 

B z d z u, 

-£l a qb x + B z d z v', 
B z d z w', 



CB z d z b' x , (8) 
(9) 

CB z d z b z ,{10) 

(11) 
(12) 

(13) 
(14) 



i'w' and p' denote perturbations in the ra- 
dial, azimuthal vertical velocities and pressures, respec- 
tively, while b' x y z represent perturbations in the magnetic 
field components. It should be recalled that q(x) is unspeci- 
fied at this stage. However the existence of a steady state is 
assumed for a given reasonable q{x) profile. Normal mode 
perturbations are assumed which are periodic in the ver- 
tical direction. A perturbation variable f'{x,z,t) is thus 
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expressed via the ansatz f'(x) exp(ert + ikz) where a is the 
temporal response while k is the vertical wavenumber. The 
incompressibility condition (1111) means that a streamfunc- 
tion tp' may be denned such that, 



u' = d z tp', w' 

Similarly the expression (fT"2]) and (fT4"]) allow for the defini- 
tion of a flux function $ in which the radial and vertical 
components of the perturbed magnetic field are expressed 
via b x = d z & and b z = —d x &. 

Utilizing these definitions, the linearized perturbation 
equations are reduced to a single one for the strcamfunction 
(wherein and henceforth the primes are dropped), 



d 2 x ip - K 2 i) = 0, 



in which 



1 



u 2 a 2 



CBik 



2 ;„2 



(15) 



(16) 



where uj az is the Alfven frequency and the spatially depen- 
dent epicyclic frequency is denned by oj(x) 2 = 2(2 — q)^Q- 
For later usage k is rewritten in the alternate form, 



= fc 2 1 



(a 2 



,2 \2 



(17) 



The solutions to this system is governed by several parame- 
ters including the vertical wavenumber k, the Cowling num- 
ber, C, the radial size of the domain (appearing below) L 
and the form and amplitude of the deviation shear profile 
Q. Taken as a representation of a small disk section it is 
assumed that there exists a minimum vertical wavenumber 
fco for disturbances since it is not reasonable to consider 
fundamental disturbances whose vertical extent is much 
larger than the disk height itself (Curry et al. 1994; and 
see Discussion). 

Before proceeding to these results it is important to 
keep in mind a number of matters. First, Knobloch (1992) 
showed that for vertical field configurations like that con- 
sidered here, together with boundary conditions in which 
individual disturbances or their combinations are zero on 
the boundaries, the character of the normal mode response 
of the MRI is that a is either real (growing/decaying expo- 
nential modes) or imaginary (oscillating modes) , but never 
complexQ A version of the argument leading to this con- 
clusion is given in Appendix [21 Thus when the character of 
the fluid response is analyzed, a 2 will be taken as e Reals. 
The relevance of this observation is that as a function of the 
system's parameters, a mode becomes unstable by passing 
through a 2 = 0, sometimes known as exchange of stabilities 
(Chandrasekhar 1961). As such, in some instances the sta- 
bility analysis of this system will focus on identifying which 
mode becomes marginal (i.e. a 2 = 0) and under which pa- 
rameter conditions. 

Secondly, the very same integral argument developed in 
Appendix [A] also leads to a general statement about the re- 
lationship between the temporal response and the structure 
of the modes in the domain, 



1 + fc 2 H (a 2 



4ft 2 , 



= 0. 



(18) 



4 It was also shown in that study that complex values of a 
can only occur if the field configuration is helical. In this case 
the instability is of a traveling wave. 



/3 is like an average radial wavenumber of the disturbance, 
2 _ Iv \9xi/}\ 2 dx 



J v \ip\ 2 dx 



(19) 



and, evidentally, (3 2 is always greater than zero, q represents 
the mode weighted average shear over the domain 



9o 



J v QH\ 2 dx 
f v \ip\ 2 dx 



(20) 



One immediate consequence of this is that if q = (i.e. 
the mode weighted average shear over the domain is zero) , 
then there is no possibility of instability. This follows from 
evaluating the criterion for marginality in (|18[) . setting a 2 
to zero and revealing 



2fllq 



0. 



(21) 



The above expression cannot be satisfied if q is zero. In 
general this relationship is not useful for direct calculations 
as the quantities therein depend implicitly upon the eigen- 
value er. However, it will be called upon to aid in inter- 
preting the results that are obtained from specific direct 
calculations in Section 6. 

5. Results 

5.1. Classical MRI modes and Channel Solutions: A Review 

The classical system is recovered for Q = 0. This means 
that k is independent of x and 



= k 2 1 



(a 2 



2n 2 q 
o 2 + ujI 



(22) 



Solutions of (|15|) are sought on a periodic domain L and 
are given in general by 



tp = Acosh(Ko.T + ■ 



(23) 



where A and <fi are constants set by the boundary conditions 
of the system. The periodicity condition that ip(x + L/2) = 
ip(x — L/2) then means that 



1 — cosh [Lk ] = 0. 



(24) 



which has the solution kqL — 2nin where n is any integer 
including zero. This condition recovers the incompressible 
limit of the classical dispersion relationship of the ideal case 
(e.g. Acheson & Hide 1973; Balbus & Hawley 1991) where 
specifically, 

{pl + fc 2 ) (a 2 + u? az ) 2 + Lo 2 k 2 a 2 - 2n 2 qokW az = 0. (25) 

in which Wq = 2(2 — qo)^! 2 , is the epicyclic frequency and 
where the radial wavenumber f3 n is 



2^2 



-; 2 



4n 2 7r 



The solutions for the temporal response, written in terms 
of respected pairs, is given by 



a 2 = a 2 Jfil) = - I . 

fc 2 



fc 2 



131 + fc 2 2 



± 



Li.'," 



0J„ 



Pl + k 2 2 



2Qlq k 2 



1/2 

(26) 
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The solutions associated with the '-' branch are termed the 
hydrodynamic inertial (HI) modes while those associated 
with the '+' are the hydromagnetic inertial (HMI) modes 
(viz. Acheson & Hide 1973). For a given mode n, the mag- 
netorotational instability occurs for the latter of these pro- 
vided that cjn > and 



2fllq k 2 



< 0. 



(27) 



The background field must be sufficiently weak for the HMI- 
modes to be unstable. 

For given values of oj 2 z the absolute minimum criterion 
required for instability to set in (i.e. for the mere possibility 
of o~ > 0) is when n — 0, i.e. kq — 0, and the condition in 
(|2"7|) predicts this to happen if, 



2Q. 2 q 



> 0, 



(28) 



(Velikhov 1959; Acheson & Hide 1973). This particular spe- 
cial mode, in which its radial structure is uniform, is gener- 
ally referred to in the literature as the channel mode which 
describes purely 2D flow with no vertical velocity. The tem- 
poral response in this case is given by a 2 n = a~ in which, 



ensure that the resulting perturbation solutions remain pe- 
riodic on the length scale L. For boundary conditions other 
than periodic ones (Sections 5.2.1 and 5.2.2), the multiple- 
time scale analysis employed is a generic procedure involv- 
ing perturbation series expansions and imposition of solv- 
ability conditions at successive perturbation orders (Bender 
& Orszag 1999). The role of the solvability conditions is to 
make sure that higher order perturbation solutions satisfy 
the boundary conditions of the system. 

The shear with zero average on the domain of length L 
may be decomposed into the Fourier series expansion 



( 2wk 

Q = eQ x ; Qi = g» sin ( — —x 



(31) 



where the parameter e <C 1 measures the overall severity of 
the shear. q n measures the Fourier amplitude of the shear 
component n and is treated hereafter as a tunable parame- 
ter. Reference to © indicates that this form for Q ensures 
that the deviation steady velocity U and, hence, the radial 
pressure gradients are zero at the boundaries x = ±L/2. 
The governing equation, with Q as given above, is expressed 
as 



2 , -0 



± 



1/2 



(29) 



When unstable (for the HMI-mode branch) this radially 
uniform mode plays the central role in the development 
of MRI induced turbulence in numerical experiments. If 
2 — > 0, then the expression for a 2 is 



cr„„ + u„ 



-fij»(2 -q ) 



±ng(2 - q ) 



-,1/2 



n 2 (2-q ) ^(2 -go) 2 



(30) 



The combination expression a 2 



(which shall appear 



on a number of occasions in the following discussion) is 
always positive for HMI modes and is always negative for 
the HM modes since the expression inside the square-root 
operation is always positive. 

It should be noted that for Q = there are no localized 
normal modes possible for a domain which is infinite in 
extent. Put in another way, there is no solution of (|15D for 
modes on an infinite domain with both (a) all quantities 
going to zero as x —¥ ±oo and (b) k is constant as given in 

(22). 

5.2. Weak shear variations on a finite domain: < Q <C 1 

In this section weak shear variations are considered whose 
average are zero on the domain L. Before proceeding a 
caveat ought to be stated. Despite the spatially periodic 
nature of the shear Q, the analysis necessary to determine 
the temporal response uses a multiple time-scale analysis. 
When periodic boundary conditions are imposed, as they 
are in Section 5.2.3, the multiple time-scale analysis is akin 
to a Floquet analysis. In this instance a long-time scale 
behavior (in the form of correspondingly weak corrections 
to the temporal response a) must be invoked in order to 



q2 / 2 i 2Q,q)z 2 Qi 

rr^ -4- i, )* 



(32) 



Solutions of this are developed in a singular perturbation 
series expansion in powers of e. Thus the ansatz, 



a 2 = a 2 + ea\ + e 2 cr 2 + • • 
«n = <L + en 2 + e 2 n 2 + 



(33) 



Given the definition of Kq given in ((22]) it follows that k 2 
K o(°o) an d 



da 2 
1 fd 2 K 2 



2 V da 4 



wherein 



Ok 2 



do 2 



o 2 +u 2 az V o 2 +uj 2 az ' (a 2 +co 2 az ) 2 



(34) 



(35) 



The possibility that a\ — must be allowed (see below) 
This is why the next order term is also included in the 
expansion. With these expansions in mind (1151) is written 
out for each power of e. To lowest order 



and to order e 



2 2^q/c (i^i 

C 0l ~ n 2 T 2 



(36) 



(37) 



The order e 2 form is relegated until later. Solutions to (f36l 
1571) are developed in the following subsections for various 
boundary conditions and restrictions at x = ±L/2. 
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5.2.1. Channel Conditions 

Channel conditions require that the perturbed radial ve- 
locity is zero at the boundaries x — ±L/2. This, in turn, 
amounts to requiring that iktp = at these positions, or in 
terms of the perturbation variables 

V'o = — 0, at x = ±L/2. 
The lowest order solution is given by 



An 



sinj m (x+L/2): 7 m = — , (38) 



CD 
"D 

=5 



3 
> 

CD 

Q 
"O 

CC 
O 
C/3 




and where Kq = — 7„. The general solution for ipo is writ- 
ten out as a sum of the linearly independent solutions ipom- 
The index m ought not be confused with the index n which 
refers to the Fourier component of the shear profile Q ap- 
pearing in (j3T1) . Stability, on the other hand, is evaluated 
by examining each individual eigenmode m. Thus all eigen- 
values a (and their perturbative corrections) will be identi- 
fied according to which eigenmode under examination, e.g. 
comi a im etc- Note also that here the index m does not in- 
clude m — as this would correspond to the trivial state 
ip = 0. 

Given the functional form of Kg found in (|2"2"|) . these 
relationships imply that <7q is equal to the functional form 

al n found in (J2B]) except with j3 n replaced by j m , i.e. ctq = kL is very large. In order for Oq X to be nearly zero it must 



Radial Position: x 



Fig. 2. The n = 1 component of the shear profile ([31]) . 
The solid line q\ = 1 and the dashed line q\ = —1. For 
problems in which channel wall conditions are imposed the 
solid line corresponds to a stable deviation profile while the 
dashed line corresponds to an unstable profile. The situa- 
tion is reversed when fixed-pressure boundary conditions 
are imposed. 



On luuuu 111 " jsim vvit-ii j~/ n L^yia^\j^x uy frm i.e. '-'q 

a 0n(ln)- The stability of these modes is now a function of 
the mode number m and, henceforth, the summation sign in 
(|3"5|) will be dropped and it shall be henceforth understood 
that when reference is made to ipo, a particular eigenmode 
mode m will be analyzed. At the next order it can be seen 
that in order to have solutions that satisfy the boundary 
condition that iktpi — at the two endpoints, the terms 
on the RHS of (|37D must satisfy a particular solvability 
condition which will relate the correction a± to the variation 
Qi. This is seen by multiplying (|3"Tf by tp given in (|3"5|) 
and integrating the result across the domain from —L/2 to 
In a general sense, then, one finds, 



be that 



2%q , cf. 



It follows that the correction 



for this mode, on, is given approximately by the solution 
to 



'ii 



^QgQ 

2(2 + 9o ) ^(2 + q Q ) 



<7i" 



(42) 



9kq 
da 2 



2 



< QlV'Om > 



°"0m + U az < n m > 



in which the bracket notation is 
o dx. 



<o>= 

< -L/2 

It follows that 

°im = -<?™(-l) m T(m), 



(39) 



(40) 



(41) 



where the last equality is established because the marginal 
mode, (Jq! ~ is under examination. The introduction of 
this profile pattern stabilizes this otherwise marginal mode 
if c/i < and destabilizes it when qi > 0. The corresponding 
shear profile for this example is depicted in Figure O It can 
be seen that the modified shear is stabilizing when the shear 
is positive near the boundaries and negative in the interior. 



5.2.2. Pressure Conditions 

An analysis of (|10p and (fT4"]) shows that the total magnetic 
plus mechanical pressure perturbations are proportional to 
the radial gradient of ip i.e., 



p/po + B z Cb z 



(43) 



in which 



T(m) = 



fig(tr, 



,2^2 



2f>2(2 + 9 ok 



Fixing the total pressure at the boundary in a Lagrangian 
sense (Curry et al. 1994) is equivalent to setting to zero the 
above expression at x — ±L/2 which, in the present case, 
means setting d x ip = at those locations. The solution 
-00m is then given by 



The denominator term of T(m) never crosses zero for phys- 
ically relevant values of the parameters so that T > 0. 

For illustration consider the m = 1 mode (which is also 
the most unstable one) near marginality in the limit where 



Vom = A m cosj m {x + L/2); 7 m = — , 



5 This procedure of applying a solvability condition actu- 
ally requires multiplying across by the adjoint solution of ipo. 
However since the eigenfunctions comprising and the linear 
operator of (|36|l are real and Hermitian, the adjoint is 



(44) 



which is very similar in form to (|38l) except that the solu- 
tion in the domain is composed of even functions and the 
series includes the mode with m — 0, which is the classical 



6 Note that this is valid since by construction the radial gra- 
dient of the steady state pressure field is zero at the end points. 
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channel mode. The correction to the temporal response for 
all values of m ^ is 



However, given the relationship between Q and Q\ it follows 
after an integration by parts that in fact, 



g m (-l) m T(m), 



(45) 



which is the negative of the solution found in (|41l) . As 
Figure [5] indicates, this result says that shear profiles in 
a channel configuration that promote stability (instability) 
play a destabilizing (stabilizing) role for configurations with 
fixed-pressure boundary conditions . 

However the m = mode is a special case which requires 
a separate analysis. In its simplest guise, the channel mode 
ipoo is a constant with respect to x. Because of the period- 
icity of Q\ it means that < QiV>oo >= and the analysis 
advanced up until this point to evaluate the correction a± 
shows that it will be zero at this order. Thus given this 
and the fact that K 2 g — the solution must be expanded 
according to, 



2 



?Aoo + #io + e 020 + 



a, 



oo 



e a 



20 



2 / ® K \ 2 



' oo 



fd ! 4\ = _ 



+ 



y 00 



(a 2 +uj 2 Y 



(46) 



> (47) 



9*4 
da 2 



'20 



2 k 



y 00 



< (d x Q) 2 > 



(53) 



Given that for HMI modes 



da 2 



>0, 



the above solvability condition says something remarkable: 
that the channel mode is always destabilized no matter 
what Qi happens to be. By contrast, values of q n may be 
chosen so that each eigenmode of the system with non- 
trivial radial structure (i.e. m > 0) may be stably influ- 
enced by the shear Q\, yet, there is no such configura- 
tion of the disturbed shear flow that can act to stabilize 
a channel mode. This result, above all others, emphasizes 
the important uniqueness of these modes. Also worth not- 
ing is that while the corrections to the growth rates of all 
modes (except m — 0) is proportional to e, the correspond- 
ing growth rate correction for the channel mode scales as e 2 : 
which means that though the channel mode is consistently 
destabilized by any periodic shear profile, its influence is 
markedly weaker. 



where a 2 the response of the channel mode is given in 
The equations for the solution of tp are in successive order 
given by 



9*000 = 0, 

9x^10 = - 



2Ci 2 k 2 Q 1 



00 • 



9^20 



{da*), 



020^00 



2nlk 2 Q 1 



(48) 
(49) 

(50) 



The first of these says that the lowest order solution which 



satisfies the boundary condition is a constant V'oo 
The next order down the solution is 

V^io = 5 " 5 QAqq, 



= A 



oo- 



at + ui 



(51) 



where 



n=l v ' 



. {2nir 
sm — —x 



\ L 
2nix 



x(-l) n cos9 n 



(52) 



In other words the solution Q is such that d 2 Q = Qi which 
satisfies the boundary condition 9 x QU=±i/2 = 0. In order 
to generate a bounded solution for ^20 a solvability condi- 
tion must be applied to the RHS of (f50|) . which is obtained 
by re-expressing tpiQ in terms of its solution as given in (|51[) 
and then integrating the equation from across the domain 
and setting the result to zero: 



da 2 



'20 



< QiQ > A. 



no- 



5.2.3. Periodic Conditions 

Solutions of tpom that are periodic on the domain L are 
given by 



-00m = A m sin j3 m (x + L/2) + (/).„ 



(54) 



where j3 m = 2mir/L and (f> m is an arbitrary phase. The 
temporal response is given by ao m (/3 m ), as defined in ([26|), 
As in the previous section, if attention is first given to those 
modes where m ^ 0, then one finds that 



2 



■ f n 
3m — 

V2 



2(f> m T(m). 



(55) 



As in the previous section, the correction a\ m for the chan- 
nel mode (m = 0) is zero and the first non-trivial temporal 
response comes at order e 2 . The procedure is exactly the 
same as found in Section 5.2.2 and its correction a 2 is 



\da 2 
with 



Q 



'20 



2 k 



< (d x Qf >, 



n=l 



(— 

\2mr 



sm 



2nn 



TT 

X+ 2 



(56) 



(57) 



Imposing periodic boundary conditions also implies that 
there are values of the coefficients of the perturbed shear Qi 
which can be chosen to promote stability for all the modes 
of the system except for the channel mode. For example, 
for the m = 1 mode and given its relative phase <pi one can 
choose q2 so that the RHS of (|55|) is negative. This can be 
done for all other modes of the system as well except for 
the m = channel mode, which is destabilized no matter 
what choice is made for the parameters of Qi, i.e. the set 
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5.3. Single shear defect on an infinite domain: localized 
disturbances 

A profile with a single shear defect located at x = is 
represented by the delta-function profile 



Q = Q Q LS(x). 



(58) 



The corresponding departures from the mean Keplerian 
shear profile U is given by 



U 



-Q n L/2, x < 0; 
Q n L/2, x>0. 



(59) 



If ip± represents the solution of (|15|) respectively for x > 
and x < 0, then solutions which decay as x — > ±oo are 



(60) 



where kq is as given in (|22[) and provided that Kq > 0. If the 
latter condition is not met, then there are no localized nor- 
mal modes allowed. The coefficients of the streamfunctions 
in both regions are equal in order that they be continuous 
across x = 0, thus A + = A- = A. Integrating (fT5")) in an in- 
finitesimal region around x = 0, and given the defect (|58[) . 
imposes a jump condition for the derivatives of ip approach- 
ing either side of the defect, 



2 k 



-QqLA. 



(61) 



Putting in the solutions for the respective domains gives 
the quantization criterion 



«Q 



• 2 +a;2 



-Q L. 



(62) 



Solutions of the above equation are constrained to two pos- 
sibilities: (i) if Qo < 0, then a 2 + ui 2 z < for a normal 
mode solution to exist while (ii) if Qo > 0, then it must be 
that a 2 + uj 2 z > for there to be a normal mode. The first 
of these cases says that normal modes are allowed only for 
those values of a 2 which are negative which are associated 
with the HI modes. The second of these cases are associ- 
ated with the HMI modes of which the MRI is a possibility. 
Thus it is an interesting result that the MRI is a permitted 
localized normal mode disturbance only if the shear defect 
is positive. If the shear defect is negative (so that there is 
an effective reduction of the shear inside the region), then 
there are no normal modes predicted. 

The solution of (|6"2")l is obtained by first taking the 
square of both sides of (JHU) yielding, 



2 T 2 



i° 2 +u 2 az ) 2 a^+ul z (a 2 +uj 2 az y 



:QoL 



(63) 



However one must be careful in interpreting the solutions 
of this simplified equation. Naively solving (1531 will yield 
twice as many solutions than are allowed. This overcounting 
is corrected by requiring that only those solutions of (|63p 
that satisfy 

Q /(a 2 +^ 2 )>0, 

are permitted or else (fB"2"|) cannot be satisfied. Solutions of 
([62| are, 



Because the term inside the radical sign is always greater 
than zero for qo > 0, the "+" branch of the above expression 
is always greater than zero while the "-" branch is always 
less than zero. Thus the "+" branch is the permitted solu- 
tion for Qo > while the "-" branch is the allowed one for 
Qo < 0. Incorporating this formally yields the expression 



-uj 2 + sgn(Q )y / ^ + 4(Q2 fc 2 L2r! 4 + ^ 



2 ) 



Rewriting the terms appearing in the interior of the radical 
sign leads to finally, 



= - \o>i. + -£) +sgn(Q ) 



-u4 {u 2 lz -2n 2 qo)+nt{kL) 2 Q 2 



1/2 



(64) 



which is similar in content to the classical channel mode 
dispersion relation (|2"9")) except for the fact that at any one 
time there exists only one branch of modes as a possible 
normal mode solution. The MRI mode only manifests as a 
normal mode when Qo > an d its growth rate is enhanced 
by the defect. 



5.4. Symmetric shear step on an infinite domain: localized 
disturbances 

A finite version of the shear profile evaluated in the previous 
section is, 



x < 



Q 



X < ~~~, 2 ) 

X > — • 

2 ' 



(65) 



This profile for Q has a top-hat structure and this de- 
scriptor will be used interchangeably with the expression 
symmetric shear step. Evidentally the value of the inte- 
gral, J_ Qdx — QoL, is the same for both forms of Q 
given in (|S"B]) and (|6"5)) . Solutions to (Ti"5l) must be developed 
separately in the three regions and appropriately matched 
across the boundaries separating the regions. Because Q 
is bounded it will suffice to match both the streamfunc- 
tions and their first derivatives across the domain. Thus, 
for x < -L/2 and x > L/2 (JT5) is 

d 2 x ^ - k 2 ^ = 0. (66) 
The solutions in the region where x < —L/2 is given by 

i>_ = A_e lKMx+L/2 \ (67) 
while for x > L/2 

1p_ = ^ +e -|«o|(x-L/2) ) ( 68 ) 

together with the constraint that Kg > 0. In the region 
-L/2 < x < L/2 ([IgD is 
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which has the two linearly independent solutions Q 

. , . . „ sinhkxl 
ipo = A cosh kx + -t>o ■ 

K 

in which 

2n 2 Q k 2 



K («o) = K Q 



SW(/!o,wo)' 
The functional form 

in which 
E (±) 



(70) 



(71) 



(72) 



<4) 



1 



(l-f 



I6f^w 



2, ,2 



O^az 





1/2" 


HI 







.(73) 



is introduced in order to make the following discussion more 
transparent. The ± designation references the HMI or HI 
modes respectively. It is also worth noting that given the 
constraint that a 2 must be real, that there is an absolute 
maximum value allowed for Kq and this is dictated by the 
requirement that the terms subject to the square root op- 
erator must remain greater than zero or else the value of 
a will turn out complex. In other words, given the reality 
of a 2 one does not expect solutions of Kq to exceed K r 
where 



2 

max 



= k 2 



i 



(74) 



In all of the numerical solutions calculated below it is 
found that, indeed, all the solutions determined satisfy 



k 2 < K max . 



Imposing the conditions that the streamfunction and 
its first derivative match across the two boundaries x 
■ 1. 2 reveals that normal mode solutions exist provided 
the following condition is satisfied 



tanh 



K 

l«ol 



coth 



K 

l«o| 



= 0. 



(75) 



The terms in the first bracket represent a quantization con- 
dition for odd-parity modes, which is to say that in the 
dimpled region ipo ~ sinh[Kx]/K, while the second repre- 
sents even-parity modes where similarly ipo ~ cos1i[kx] in 
the same region. 

Attention is first given to the even-parity mode in which 
the quantization condition, 



coth 



L 



K 

l«o| 



must be satisfied for some value of kq. If k were real and 
greater than zero, then the LHS of the above relationship is 
always positive. Since Kq is real and positive, if the modes 
are to be localized, then the RHS of this expression is always 
less than zero. In this case it means that no solution exists 
for k real. On the other hand solutions do exist if k 2 < 0. 



7 The second of the solutions appearing in (|70p is written in 
the form displayed in order to make sure that both linearly in- 
dependent solutions are represented even in the event K = 
(Friedman 1956). 



Rewriting k — ik then the quantization condition above 
becomes 



cot[KL/2] = 



l K o| 



(76) 



now with the constraint for the existence of localized normal 
modes is re-expressed as 



K 2 



(«o) 



2n 2 Q k 2 

S( ± )(/«o,wo) 



> 0. 



(77) 



The strategy for obtaining a solution is as follows: find the 
values of kq which solve (1761) , subject to the constraint given 
in (|77[) . and then given the solution kq find a 2 through 
the relationship ((72|) . The same procedure will be imple- 
mented for the even-parity solution below as well. Solutions 
of the quantization condition will be sought graphically in 
the range of values for kq € (0,K maa: ). Inspection shows 
that the constraint (fTTf is violated if Qo < and £w > 0, 
or if Q Q > and < 0. Since EW > it means that 

there are no HMI normal modes as allowable solutions if 
Qo < 0. HI normal modes are associated with the E^ - ) 
branch of solutions and since this is less than zero except 
for very special circumstances, one can generally expect HI 
normal modes for Qo < and for them to be stable. The 
aforementioned special circumstance occurs if there are so- 
lutions to (1751 in which kq > k since, in that case, > 
and HI modes may exist as well for Qo > 0. 

The even-parity mode limits to the single-defect con- 
figuration discussed in Section 15.31 for small values of the 
parameter L. Analysis of (1761) in this limit shows that, 



2 

kL 



\ K o\ 



which, after restoring the definition of k and some rearrang- 
ing, becomes 



«o = TT 



L /2^Q fc 2 
2 \V 2 



•us 



where the definition of has been restored. For values of 



Qo such that fioQo > {o 2 - 
reduces to 



z )n^/k 2 the above expression 



«o = 



n 2 Q Lk 2 



which is the quantization condition for the single-defect 
problem on the infinite domain, (f6"2"j) . The range of val- 
ues for Qo and L for which this limiting form agrees with 
the actual result of the quantization condition is shown in 
Figure |31 where u> 2 z is set to the value 2Hq qo in order to 
compare the temporal response against conditions in which 
the channel mode is marginal in the classical theory, see 
(|28[). The main result is that the temporal response of the 
defect profile of the previous section closely resembles the 
temporal response of this step profile here in the limit where 
kL is small less than one. The limiting form represented by 
the defect profile becomes a poor representation of the top- 
hat profile when kL exceeds order 1 values, however, the 
agreement tends to be much better for the HI. 

Additionally, for small values of kL there is always at 
least one normal mode expected, whether it be the HI or 
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HMI mode. However, as the horizontal length of the sym- 
metric step profile increases, the number of permitted nor- 
mal modes increases as well (Figure 2]). Similar behavior 
is predicted for modes which are localized in the vertical 
direction of the disk (Liverts & Mond 2009). As the num- 
ber of permitted normal modes increases their ko values 
will cluster (countably as kL — > oo) around the maximum 
value 



(«§) 



Qo 



12,. ,2 



7,8) 



Odd-parity modes are interesting since they have no 
analogues in the single defect profile studied in the previous 
section. Examining the quantization condition for it, 



tanh k 



L 



K 



shows that there are no real solutions of k that satisfy the 
constraint. Utilizing the expression k = ik, as used previ- 
ously, the quantization condition takes the form 



tan[fcZ/2] 



K 

Kl 



(79) 



Inspection of the condition ([79)1 shows that for there to 
exist a normal mode k must be greater than n/L as the 
functional form ta,n[kL/2]/k is positive for values of k G 
(0, tt/L). This places a constraint upon the minimum value 
of L that permits this mode to exist and this obtained by 
solving 



7T 

kL' 



(80) 



and, using the definition for k found in f|71[) one finds that 
the minimum value for L must satisfy, 



9 2 

4Q 



sgn(Qo 




1/2 



(81) 



Note that the inclusion of the expression sgn(Qo) reflects 
how HMI/HI modes associate with the sign of Qo so that 
the expression appearing above encapsulates both possibil- 
ities. Odd-parity modes bifurcate into existence for L — >• oo 
as Qo grows from zero and therefore it means that these 
modes are not expected to play a role for values of the de- 
fect which are small as this would entail very large values 
of the top-hat region. The temporal response of these odd- 
parity modes are otherwise similar in most every respect 
to their even-parity counterparts and, as such, no further 
analysis of their properties is pursued here. 

It should be noted that as in the previous section un- 
stable localized MRI normal modes (HMI-modes) can exist 
only if the step function in the shear is greater than zero. If 
the step is negative (indicating a region of weakened shear), 
then there are no HMI types of localized normal modes ad- 
mitted by the configuration and, hence, there is no possi- 
bility of normal mode unstable disturbances. 



10 J 



10 u 



10" 



(a) Temporal Responce, = 0.6 (HMI Mode) 



10" 



10 u 



10' 



kL 



(b) Temporal Responce, = -0.6 (HI Mode) 



-10" 



-10 



0.87 



0.5 



1.5 
kL 



2.5 



Fig. 3. The temporal response resulting from the quanti- 
zation condition (I76D as a function of kL (solid line) com- 
pared to the quantization condition (1521 for defect profile 
(dashed line). For both plots uj 2 z = 2Sl§g with Qo = 3/2. 
(a) For Qo > only HMI normal modes are permitted. The 
agreement between the two profiles severely breaks down 
for kL exceeding 2. (b) For Qo < only HI normal modes 
are allowed. The agreement between defect theory and step 
profile is better for HI modes. 



5.5. Single shear defect on a periodic domain 

Next a shear defect of the form 
Q = -Qq + QqLS{x), 



(82) 



is considered on the periodic domain —L/2 < x < L/2. 
This form is very similar to the one considered in Section 
I5.3l except the integral of the shear over the periodic domain 
is zero, i.e. < Q >= 0. The effects of shear defects with zero 
mean may be examined. The governing equation for ip is 



2VL 2 k 2 Q L8{x) 



where k 00 is given by 



= fc 2 1 



2n 2 (q - Qo) 



(83) 



(84) 



Similar to what was done in the previous sections, one may 
express a in terms of k 00 via the expression and, for 
the sake of completeness, this is written out explicitly, 



-L = s (±) (- 
i 



2 \ _ 
oo» W ooJ = 



2 (l - lk) 



UJ 2 ± 



no 
k 2 



-,1/2" 

(85) 



in which the shear modified epicyclic frequency is defined as 
co 2 = 20§(2 — <7o + Qo)- Solutions to (1551 must be developed 
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Fig. 4. Plot showing the increase of allowed normal modes 
with domain size L. In all plots lu 2 z — 2J7q<7o, <Zo = 3/2 and 
Qo = 0.6 (HMI mode permitted). The dashed line is k/|ko| 
and solid line is cotkL/2. Panel (a): kL = 2tt corresponds 
to one normal mode at kq = 0.50. Panel (b): kL = Yin 
supporting four normal modes kq = 0.32,0.465,0.54,0.56. 
All modes correspond to values of a 2 > 0. « (max) = 0.576. 



separately for either side of x = and matched. Solutions 
for the streamfunction in which ip(—L/2) = ip(L/2) and 
d x ip\-L/2 = dxii\L/2 (i- e - that they are periodic) are given 
by 



ip ± = A cosh 



XT- 



L 



(86) 



where ip ± is the solution for the region x € (0, ±L/2). 
Continuity of the stream function is ensured by construc- 
tion of the form of the solution in (|56"|) . The streamfunction 
must show a jump in its first derivative following the same 
arguments found in Section 15.31 Integrating (1831) in a small 
region around x = results in the same condition for the 
jump in d x ip found in (1611) . Putting in the forms for ip ± 
found in (|5o) leads to the quantization condition 



tanh 



n 2 Q k 2 Q L 

= {±) Ko, <) 



(87) 



It should be kept in mind that the quantities as 
appearing on the RHS of (1571) are different functions of 
k 00 depending upon whether one is considering HMI/HI 
(±) modes, respectively. Thus when solutions of the quan- 
tization condition are sought, separate attention will be 
given depending upon which type of mode is of interest. 
Furthermore, consideration of the solutions to ((87)) must 
be done by restricting attention to cases where kqq is either 
real or imaginary. Additionally from previous arguments 
showing that a 2 must be real, there will be a maximum 
real value possible for kqq 



2 12 

Krmn^. \LO 



= k 2 



Values of k 00 which are greater that K max (lu 2 ) would lead 
to complex values of a 2 as an inspection of (l8"B"j) shows that 
if kqq does exceed this maximum value, then the term inside 



the square-root operator becomes negative which would 
lead to a complex value of a 2 . This restriction upon the 
value of k qo helps in finding solutions for it. 

The limit where Qo — > recovers the classical limits 
discussed in Section 15.11 In that case channel modes, i.e. 
those modes for which kq = 0, are the most unstable. In 
the current setting this corresponds to n m — when Qo 
— 0. Additionally, for uj 2 z = 2S!q(7o the channel mode is 
exactly marginal for Qo — 0, viz. the arguments leading 
to (f2"5| . This value for ui 2 z will be assumed so that a con- 
trolled analysis may be carried out in which one can track 
how varying Qo affects the stability characteristics of an 
otherwise marginal channel mode. 

Figures [5][6] show the graphical solutions of ((87|) with 
some details of the results. The main result is that the mode 
that is identified with the classical channel mode in the 
ideal limit always exhibits exponential temporal growth no 
matter what amplitude the shear defect takes. Figure [7] 
shows that the temporal response of the classical channel 
mode is maximal for Qo = and that as Qo moves away 
from zero instability is predicted. The mode exhibits an 
exponentially decaying spatial character (k 2 > 0) for Qo > 
0, while it is weakly oscillatory (k 2 q < 0) for Qo < 0. This 
general trend persists for a wide range of kL values. As 
such, the results depicted in these figures can be taken to be 
qualitatively representative. The figures also indicate that 
the unstable HMI mode can also turn into an unstable HI 
mode if the amplitude of Qo becomes large enough. 

5.6. Single Shear Defect on Finite Domain: Channel 
Boundary Conditions 

This section is concerned with the response of disturbances 
to a Q configuration that is exactly like that considered in 
Section 15.51 where the difference now is that the tj) distur- 
bances are set to zero at the boundaries x = ±L/2. This is 
equivalent to enforcing zero normal velocity at the bound- 
aries. Thence, ip in this case is 



sinh [k 00 (it |)j 
sinh k 00 L/2 



(89) 



where just as in the previous section, ip ± is the solution for 
the region x € (0, ±L/2). Evaluating the jump condition 
at x = leads to an analogous quantization condition like 
that in (|57|) . namely, 



k 00 coth 



K °° 2 



fl 2 k 2 Q a L 



n 2 k 2 Q L 



£(±) 



(90) 



in which £w (k oo; u 2 n ) and k 00 and lo 2 b are exactly as they 
are found in Section [5751 It can be seen here, as before, that 
there there are never real solutions of k 00 for HMI modes 
if Qo < and, likewise, no such solutions for HI modes if 
Qo > 0. 

In Figure [5] the stabilizing behavior of the shear defect 
on the HMI modes is depicted for the most unstable mode 
of the system. The results from the weak shear limit dis- 
cussed in Section l5.2.1l are used as a reference point to help 
interpret the following. From the classical theory discussed 
in that section it was established that the most unstable 
normal mode allowed is the one in which m = 1, in turn 
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Fig. 5. Graphical solutions for the quantization condition 



([87]). Attention restricted to real values of k 00 marking 
the mode which becomes the classical channel mode in 
the limit Qo — > 0. Each panel depicts the three functions 
k 00 tanh k 00 L/2 (thick line), and (dashed/thin lines). 

In all plots kL = n and uj 2 z = 2{l 2 t q . Panels (a) and 
(b) show that HMI modes are admitted as the solid curve 
crosses the line while for the latter two it crosses the 
line corresponding to J 7 ^-* meaning that HI modes are se- 
lected. The solution k 00 and temporal response for each: 
(a) Qo = 0.1, k 00 = 0.27, a 2 = 0.005 ,(b) Q = 0.5, k 00 = 
0.65, a 2 = 0.14, (c) Q = 2.0, k 00 = 1.33, a 2 = 3.36 and 
(d) Qo = 3.5, k 00 = 1.37, a 2 = 6.73. Notice how the chan- 
nel mode goes from being an unstable HMI mode to an 
unstable HI mode as Qo is increased. 



Fig. 6. Similar to Figure [S] except negative values of Qo 
are investigated. In all plots kL = ir and u> 2 z — 2f2o9o and 
qo = 3/2. Panel (a): the left plot shows that the channel 
mode admitted is a stable HI mode with k 00 = 0.92 and 
(T 2 = —5.30. The right plot in Panel (a) shows the first two 
overtone modes. The first HMI mode is unstable with k 00 = 
0.69i and a 2 = 0.88 while the next HI mode admitted has 
k 00 = 1.32i and a 2 = —2.88. Note that in these cases k 00 
is imaginary. Panel (b): stable HMI mode with k 00 =1.21 
and a 2 = —8.30. For the first two overtones, the HMI mode 
is unstable with k 00 = 0.78i and a 2 = 2.02 while for the HI 
mode, k 00 — 1.08z and a 2 = —2.47. For all values of Qo < 
there is at least one unstable HMI mode with imaginary 



implying 71 — tt/L. Reference to (|2~T|) shows that this mode 
is marginally stable when 

^«= 2 "fe° fc2^+V ( 91 ) 

(This is obtained by replacing /?„ with 71 in that expres- 
sion.) The response of the modes for three different values 
of kL are depicted. With the above choice for u 2 z all modes 
are marginal at Qo = 0. For the range of parameters inves- 
tigated and shown in the figure, k 00 must be imaginary in 
order for the quantization condition to be satisfied. Thus 
the values for k 00 are written as iH 00 . Inspection of the re- 
sults show that stability is promoted only for values of Qo 
which are negative. The range of values of Qo for which sta- 
bility is predicted decreases for larger domain sizes, i.e. as L 
is made larger the window in Qo for which disturbances are 
stable decreases. For those ranges in Qo in which stability 
is predicted, e.g. for Q € (—0.5,0) for kL = 2n with co 2 z 
given in (|91l) , the temporal response of the other overtones 
was checked and all were found to have a 2 < 0. 



6. Discussion 

This study is devoted to examining the axisymmetric in- 
viscid linear response of a shearing sheet environment 
threaded by a constant vertical magnetic field for an array 
of different shear profiles. The motivation for considering 
this setup are restated from the Introduction: 

1. Numerical experiments of low magnetic Prandtl number 
flows appropriate for cold disk environments are outside 
current computational reach. 

2. Current numerical experiments seem to show the trend 
that as P m is lowered so is the corresponding angular 
momentum transport. The reasons for this is not yet 
clear. 

3. Theoretical analysis of laboratory setups evaluating 
the response of the MRI subject to a vertical back- 
ground field find that transport is reduced with P m . 
Furthermore this reduction emerges through the modi- 
fication of the basic shear profile (e.g. Taylor-Couette) 
inside the experimental cavity which, in turn, is driven 
into place by the MRI itself. The final condition is a 
stable pattern state. 
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Behavior of Channel Mode with Defect Amplitude 




Q Q : Defect Amplitude 



Fig. 7. Behavior of the classical channel mode as a func- 
tion of defect amplitude Qo for kL = ir and uj 2 z = 2Q,Qqo, 
qo = 3/2, corresponding to a marginal (a 2 = 0) channel 
mode at Q = 0. Top panel depicts k 2 q as a function of 
Qo in which the classic channel mode has k 2 = 0. Bottom 
panel shows the temporal response which predicts exponen- 
tial growth/decay for all values Qo ^ 0. The channel mode 
goes from being a HMI mode to a HI mode at Qo ~ 1-87. 

4. In the final pattern state reached by the model labo- 
ratory system the modified shear profile's amplitude is 
independent of P m . This is not true of the remaining 
fluid and magnetic quantities which scale as with 
typical values A w 1/2. As P m is made very small one 
would observe a fluid state characterized by a uniform 
background magnetic field and an azimuthal flow show- 
ing deviations from a Keplerian state. 

5. The main mode of instability and driver of turbulence 
in the numerical experiments of the shearing box is the 
channel mode which is a special response with no radial 
structure. This mode precipitates a secondary instabil- 
ity which is understood to lead to a turbulent cascade. 
There is no channel mode allowed in a setup appropriate 
for a laboratory experiment. 

The saturation hypothesis is then the following: is it possi- 
ble that the fluid response in a shearing sheet environment 
characterized by very low P m stabilizes itself by modify- 
ing the background shear as suggested by the calculations 
appropriate for laboratory conditions? Of course, there are 
some constraints placed upon the way in which the shear 
may be modified in that the gravitational field predomi- 
nantly "forces" there to be a Keplerian profile (Ebrahimi 
et al. 2009). However by analogy to the situation for the 
laboratory setup, there may come into play undulations of 
the shear whose average is zero over some length scale L 
(e.g. Figure 1). 

Of course, the main difference between the laboratory 
setup and numerical calculations of the shearing box are 
their respective boundary conditions. The numerical exper- 



Behavior of Least Stable Mode 



kL = it 

kL = 2it ^s, 
kL = 4it 








kL = 2n \ 

\ kL = 4it 
kL = it >v \ 









-2 -1.5 -1 -0.5 0.5 



Q Q : Shear Defect Amplitude 

Fig. 8. The behavior of the first unstable (HMI) mode for 
the problem with channel walls (qo — 3/2). The values of 
ik 00 — k 00 and the corresponding temporal response a 2 are 
given as a function of the defect amplitude Qo for three 
values of the parameter kL. The top panel shows k 00 . The 
bottom panel shows the temporal response. The value of 
uj 2 z for each curve is set so that the temporal response is 
marginal at Qo = - details of this is described in the text. 
The value of k 00 for Qo = are shown by the V symbol in 
the top panel. 



iments adopt periodic boundary conditions while the cal- 
culations for the laboratory experiments require no normal 
flow conditions at an inner and outer radius and the former 
of these permits channel modes to come into play. 

Thus, what has been done here is a simple model calcu- 
lation in which a shearing box configuration with a constant 
background vertical magnetic field is tested for stability for 
a number of different shear profiles. The perspective taken 
is that if there exists a shear profile (say, whose average 
over some length scale is zero) that is stable to the MRI, 
then it would strongly suggest that low P m flow configura- 
tions might saturate the MRI by driving into place a new 
azimuthal flow profile which leads to saturation. 

What is found here is something very interesting. For 
flow configurations in which no-normal flow conditions are 
imposed at both radial boundaries (Sections I5.2.1I & I5.6I) 
there exists a modification to the shear in which the most 
unstable MRI mode is stabilized. However, for flows in 
which periodic boundary conditions (Sections I5.2.3I & I5.5p 
or fixed Lagrangian pressure conditions (Section I5.2.2[) are 
imposed there are no modulations of the shear which shut 
off the instability. In these cases, all modifications of the 
shear (with zero mean on some radial length scale) seem 
to enhance instability of the channel mode. For example, if 
the given channel mode is marginally stable to exponential 
growth, then introduction of the modified shear, no matter 
what is its non-zero amplitude, seems to destabilize it as 
the results of Section [5.2.3I and Section [5~51 indicate. 
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This strongly implies that in those numerical experi- 
ments in which a shear is threaded by a constant vertical 
magnetic field, the shear modification/MRI-stabilization 
process may only apply for those flow conditions in which 
there is no normal flow on either one or both radial bound- 
aries of the system. 

On the other hand, for a similar magnetized fluid con- 
figuration described by periodic boundaries which there- 
fore permits channel modes to exist, modifications of the 
shear appears to further destabilize the channel mode. If 
one were to consider a flow configuration (e.g. the local 
Keplerian profile in a shearing box threaded by a constant 
background field) that admits an unstable channel mode, 
then there seems to be no modification of the shear with 
zero mean that would saturate the growth of the chan- 
nel mode since all such modifications further enhances its 
destabilizing influence. The emergence of a pattern state 
leading to saturation like the one envisioned by the hy- 
pothesis would seem to be out of the question under these 
circumstances. Moreover this suggests that the tendency for 
angular momentum transport reduction with P m , as indi- 
cated by numerical experiments of the shearing box, is due 
to some other process or processes (which are not discussed 
here) . 

These issues deserve further investigation as the con- 
clusions reached here were done so utilizing very simple 
functional forms for the shear profile Q. Namely, the analy- 
sis performed in the previous sections either assumed weak 
values of Q, in order to use singular perturbation theory 
techniques, or that Q is described by delta-functions (see 
also below). Use of these functions, which greatly facilitate 
analysis, have offered some insights into the nature of the 
responses of these modes. Of course, further investigation 
would require considering more physically realizable profiles 
and to check the robustness of the results they yield against 
the ones offered by use of these more simplified forms. As 
such the results reached here can serve as guideposts for 
further enquiries along these lines. 

In the following some reflections are presented on other 
issues pertaining to the implications of the calculations per- 
formed here. 

6.1. On the use of delta-functions 

A comparison between the responses of a single delta- 
function shear profile on an infinite domain (Section 15 -3|) 
and a top-hat profile for Q (Section 15. 4| show that the 
results of the two configurations are qualitatively similar. 
Although the symmetric step function allows for more nor- 
mal modes to exist as its radial extent is increased - a fea- 
ture which is missed if one uses the delta-function prescrip- 
tion - the qualitative flavor of the response is captured by 
its use nevertheless. 

Another deficiency found in using the defect prescrip- 
tion is that for the top-hat profile studied in Section 15.41 
there are both even and odd parity normal modes possible 
whereas only the even-parity mode is admitted in the de- 
fect model. The odd-parity mode has no counterpart in the 
single defect model. It was also shown that the odd-parity 
mode may exist if the radial domain of the symmetric shear 
step is large enough as (15Tj) indicates that the minimum 
size of the domain required for this mode to exist has the 
proportionality L ~ 1/ \/Qq . For small values of the defect 
amplitude the minimum size of the shear step required for 



this odd-parity mode to exist is so large that one does not 
expect that this should play any physically relevant role in 
the dynamics if one is interested in the response to nar- 
row regions of altered shear. It is with some confidence to 
suppose that it is reasonable to represent the dynamical 
response of slender regions of modified shear through the 
use of delta-function defects and this is the motivation and 
justification for its use in Sections 5.5 and 5.6. 

6.2. Channel modes 

In the small Q theory of Section 15.21 it was shown that 
channel modes exist for periodic boundary conditions and 
for conditions in which the total pressure perturbations are 
zero at the two boundaries. Conversely, if either one of the 
two boundaries force the velocity perturbation to be zero 
there, then there is no channel mode permitted as a normal 
mode solution. In light of the previous results it is impor- 
tant that the boundary conditions appropriate for disks be 
properly ascertained in order to assess and better under- 
stand how activity in them are driven. 

Channel modes in the classical limit have no radial 
structure. In Section [53] the development of this mode was 
studied as a function of the defect amplitude Qo and it is 
shown that though the channel mode develops some amount 
of radial structure, as evinced by the non-zero value of its 
wavenumber k 00 , its tendency to be unstable survives and is 
even enhanced. In Regev & Umurhan (2008) it was postu- 
lated that the channel mode might disappear as a natural 
mode of the system if some amount of radial symmetry 
breaking were introduced into the governing equations of 
motion, for instance, in the form of a radially varying shear 
profile. However, the results of this section show that the 
classical channel mode indeed persists in its existence de- 
spite the introduction of a shear profile which deviates from 
the background constant shear state with zero average over 
some length scale. These conclusions are consistent with 
similar results regarding the nature of the MRI in more 
global contexts (Curry et al. 1994). 

6.3. On the existence and number of localized normal modes 

In Sections 15. 31 and T5 .41 the existence of localized modes was 
examined. Localization is understood in this case to corre- 
spond to modal disturbances which show exponential decay 
as the radial coordinate x — > ±oo. Thus, by construction, 
incoming/outgoing waves are excluded from consideration. 
The existence of localized normal modes depends on the 
relative sign of the shear defect /step profile Qo: HMI modes 
exist if the shear in the defect/top-hat region is stronger 
than the background (Qo > 0) while HI modes exist if the 
shear is relatively weakened (Qo < 0) - summarized here 



Given the mathematical structure of the ode describing ax- 
isymmetric disturbances, the radial eigenfunctions that are so- 
lutions have either strictly oscillatory or exponential radial pro- 
files. Thus the possibility of decaying oscillations is ruled out. 

9 The mode existence dependency upon the sign of Qo was 
checked for other shear profiles (for instance, by replacing the 
top-hat with a truncated parabolic profile). Aside from differ- 
ences in the magnitude of the temporal response and the mul- 
tiplicity of modes allowed, the aformentioned sign dependency 
still holds. 
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as 



Qo > < — > hydromagnetic inertial disturbances, 
Qo < i — > hydro — inertial disturbances. 

Given that the former mode type leads to the MRI, it 
means that if the shear in the defect zone is weakened, then 
there are no localized normal modes which can go unsta- 
ble as only hydro-inertial modes are supported (which hap- 
pen to be stable except for very extreme circumstances). 
Of course, as an initial value problem, it does not mean 
that there are no hydromagnetic inertial modes permitted 
but, instead, the mode probably has some type of algebraic 
temporal growth/decay associated with it which cannot be 
described by the usual normal mode approach. This needs 
to be further examined by studying the associated initial- 
value problem. Similar features are known to exist for ver- 
tically localized MRI modes in disks (Liverts & Mond, 
2009) in which the response of the initial value problem can 
exhibit exponential growth with an amplitude supporting 
some algebraic time-dependence. Growth can initially be 
very small and it can take a very long time for a given dis- 
turbance to reach the same amplitude as expected for its 
counterpart normal mode. A conjecture is that the meaning 
of the absence of normal modes in the problem considered 
here may have similar attributes to the features associated 
with the initial- value problem examined by Liverts & Mond 
(2009). 

Irrespective of this outcome it should be remembered 
that the dichotomy observed here for modes on an infi- 
nite radial domain is a consequence of the imposition of 
exponential decay of the disturbances. This distinction dis- 
appears once waves are allowed to enter and exit from the 
"infinite" boundaries. 

Lastly it is also noteworthy that for the localized prob- 
lems considered here, when normal modes are permitted 
they usually come as a finite set. This is in contrast to 
other localized flow problems considered elsewhere in the 
literature. For example, in the problem of compact rotat- 
ing magnetized jets in an infinite medium considered by 
Bodo et al. (1989) a countably infinite number of normal- 
modes are predicted for that system. This has some origin 
in the uniformity of both the rotation and magnetic field 
profiles of the jet considered in that study. By contrast, 
there are examples in geophysical flows in which the num- 
ber of normal mode disturbances are discrete and finite. A 
recent example can be found in the work of Griffiths (2008) 
where the stability of inertial waves are studied in a strati- 
fied rotating flow with strong horizontal shear. For strongly 
peaked forms of the potential vorticity only a finite set of 
normal modes are predicted for the system. Such finiteness 
of the number of normal modes permitted is not so un- 
usual in systems in which background states, like the shear 
or stratification, have strongly peaked functional forms like 
they are in the cases studied in this work. 

6.4. Stabilization in general and an interpretive tool 

The absolute minimum condition that must be met for any 
normal mode of the system to be unstable in a configura- 
tion with constant shear qo is given by the Velikhov cri- 
terion (l28|) : if the square of the Alfven frequency (wfj is 
greater than the shear by 2H,qQq, then none of the HMI 
modes of the system can be unstable. For individual modes 



in the same system with some amount of radial structure 
the criterion is given by (|27l) , where 0^ is the square of the 
radial wavenumber of the disturbance. As can be seen, for 
fixed values of ui 2 z there are two ways in which stability 
may be promoted - either by weakening the shear or by in- 
creasing the radial wavenumber. In the problem where the 
shear is constant, these quantities are set from the outset 
as parameters. 

At the end of Section 4 there is a series of arguments 
for arbitrary shear profiles which lead to the identity found 
in (|18l) . which is like a generalized eigenvalue condition. 
This relationship is the analog of the dispersion condition 
of the classical limit (uniform q) contained in (|25p F°l Care 
must be adopted before interpreting this expression as an 
eigenvalue condition since a appears implicitly in the ex- 
pressions for q and /3. Nonetheless, a parallel analysis shows 
that stability is promoted if 



2^qk 2 
fi 2 + k 2 



>0, 



(92) 



where the effective wavenumber f3 and shear q, given in 
(fj"9l20j) are, 



/5 2 



J v \d x ip\ 2 dx 
J v \ip\ 2 dx 



q = qo 



J v QM 2 dx 
J v \*p\ 2 dx 



T> is the domain and %j) is the radial eigenmode of the distur- 
bance in question. The stability condition for general q may 
be rationalized in the same way as the criterion is under- 
stood for the classical case with the replacements, f3 n — > /3 
and go ~~ ► 9- With these aforementioned observations and 
caveats in mind, one may use (|92p as an interpretive tool to 
aid in understanding the reasons for stability /instability in 
a particular problem. Thus, stability is promoted if (i) for 
fixed effective wavenumber, the effective shear is reduced 
or (ii) for fixed effective shear the effective wavenumber is 
increased, or in more general terms (iii) if the combination 
results in an overall reduction of the quotient 

<1 



/3 2 + k 2 ■ 

If one considers the results concerning the problem with 
periodic boundaries in Section 15.51 then the persistent in- 
stability expected for channel modes can be understood in 
terms of this criterion. In the classical limit where the de- 
fect amplitude is zero then (3 for the channel mode is zero. 
As the defect amplitude is raised /3 increases suggesting 
that this effect has a stabilizing influence, e.g. see Figure 
IH1 However, as the defect amplitude is raised away from 
zero the effective shear increases no matter what the sign 
of Qo and, to effect, this increase always sufficiently over- 
powers the increase in /3 2 so that the end result is further 
instability. In the same figure is plotted the quantity 



W 
Wq 



Q Pi 



qo P 2 + k 2 ' 



(93) 



which, for the generalized Velikhov criterion, is a measure 
of the ratio of the destabilizing term for the shear profile q 
compared to its nominal value when the shear is only qo: for 
zu/wq > 1 one expects enhanced destabilization against the 



10 The classical limit is obtained when the replacements (3 
p n and q — > qo are made. 
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Fig. 9. The behavior of the unstable channel mode for the 
problem with periodic boundary conditions examined in 
Section 1531 kL — n and lu 2 z — 2q with q = 3/2. Plotted 
is the temporal response, a 2 , together with its correspond- 
ing effective wavenumber and shear ft and q = qo + Q (Q 
plotted), and the stability measure for this mode w/zuq. All 
values of the shear defect Qo correspond to increasing the 
effective shear to above the background value Q(Qo) > 0. 
This increase in q always overpowers the stabilizing influ- 
ence of an increased ft 2 . Note, the values of a 2 have been 
scaled by an arbitrary factor in order to facilitate clear com- 
parison. 



Fig. 10. The behavior of the least stable HMI mode for 
the problem with channel walls examined in Section 15.61 
kL = 7r and lo 2 z = 2q /(k 2 + ft 2 ) with ft\ = tt / L and 
qo = 3/2. Plotted is the temporal response, a 2 , together 
with its corresponding effective wavenumber and shear ft 
and q — qo + Q (Q plotted), and the stability measure 
for this mode w/vjq. Stability occurs for Q £ (—1.32,0). 
In this window the effective shear is less than the back- 
ground value, (q — qo = Q < 0) for only part of this 
range, i.e. for Qo € (—0.5,0). Stabilization in the range 
Qo G (—1.32,-0.5) is achieved by a sufficient increase of 
the effective wavenumber so that the overall effect upon the 
stabilization parameter ■uj/wo is for it to be less than one. 
cr 2 has been similarly scaled as in the previous figure. 



classical mode while for vj/wq < 1 comes with enhanced 
stabilization. In Figure [9] this quantity is plotted for one 
of the channel modes discussed in 15.51 Noting that in this 
case ft n — (i.e. that its radial wavenumber is zero in the 
ideal case) the radial profile of the mode develops as the 
amplitude of the shear defect increases no matter what its 
sign. The stabilizing rise in ft is accompanied by the desta- 
bilizing increase of the effective shear, that is to say Q > 0, 
so that the quotient of their respective influences results in 
vj/wq > 1, indicating instability. 

Similarly plotted in Figure [TU] is the influence Qo has 
upon ft, q, and the stability measure vj/vjo, for the sin- 
gle defect problem on a finite domain studied in Section 
15.61 The presence of walls filters out the channel mode 
and the resulting flow potentially can support shear pro- 
files which can stabilize the least stable mode. The Figure 
shows that within the range of (negative) values of the 
defect amplitude Qo for which the mode does not expo- 
nentially grow/decay the effective shear is less than the 
background shear state only for part of the stable range. 
It means to say, then, in that range of parameter values 
in which the effective shear is greater than the background 
shear, yet there is still stability of the mode, stabilization 
is brought through an increased effective radial wavenum- 
ber which over-compensates the increased destabilization 
brought about by an increase in the effective shear. 



6.5. Final reflections and an implication 

The fundamental basis of this study is the assumption that 
in low magnetic Prandtl number flows for model environ- 
ments like that considered in this study, the only noticeable 
outcome of the development of the MRI is to alter the basic 
shear profile. All other quantities, such as magnetic fields 
and the radial and vertical velocities, have saturated pro- 
files which scale as some positive power of P m so that for 
P m <C 1 their contributions become negligible. The basis 
of this assumed trend derives from the theoretical calcula- 
tions of laboratory setups discussed in the Introduction. It 
is therefore assumed here that a similar trend may appear 
for configurations in shearing boxes with periodic radial 
boundary conditions as well. The perspective taken here 
with regards to those configurations is that if this is in- 
deed the case, then one can (in principle) test the stability 
of a wide variety of shear profiles which show deviations 
(constrained to have a zero average over some length scale) 
from the basic Keplerian profile - configurations which are 
akin to those predicted from the calculations of laboratory 
setups. If a stability calculation shows that there exists a 
shear profile which is stable to these axisymmetric distur- 
bances, then the shear-modification/mode-saturation hy- 
pothesis should become a serious candidate explanation. 
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However, the calculations done in this study indicate 
that for periodic shearing box environments there is always 
an instability of the basic channel mode no matter what 
shear profile is assumed (satisfying the aforementioned con- 
straints). Unlike the situation for laboratory setups, where 
saturation can be achieved by altering the basic shear pro- 
file, saturation (if present) in shearing box environments 
with periodic boundary conditions is likely not due to this 
mechanism. 

This further implies that the nonlinear response of sys- 
tems supporting the MRI stongly depends upon the bound- 
ary conditions employed which is not so surprising as many 
physical systems in Nature exhibit this sensitivity to their 
boundaries. Thus it seems to this author that some care 
must be taken before one equates the results of laboratory 
experiments to those physically related/ analogous systems 
for which the experiments are meant to represent (cf. Ji et 
al. 2006). This cautionary note is not to diminish the value 
of one over the other, rather, it is intended to emphasize 
that there are many subtle and, at times, conflicting fea- 
tures between such systems that must be clearly understood 
before any firm conclusions are reached. 
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Appendix A: Proof that a 2 is never complex 

(|15p is rewritten into the form 

20 2 Q 



a* - 

20g(2 - 
by supposing that 



^ = o, 



where Wq 



9o J 



The nature of a 2 will be assessed 



a + UJ az =a + ib 



where a, b G Reals. It will be shown below that b must be 
zero. 

The boundary conditions on ip are one of the following 
possibilities: (i) ip goes to zero on the domain boundary 
3D provided the latter is finite, (ii) ip and its derivative 
go to zero on the domain boundary if the boundary tends 
to infinite distances, (hi) ip is periodic on a length scale L 



in the domain V. Replacing a 



with a + ib in the 



equation for ip, followed by multiplication by the complex 
conjugate of ip, (i.e. ip*), integrating the result and applying 
the boundary conditions gives 



h + ih = 0, 
where 



(A.l) 



h = ~ [ \d x iP\ 2 dx 



2n 2 Q 



b 2 



2, ,2 



b 2 



\iP\ 2 dx{k.2) 



and 

h = bk 2 



v 



lo 2 - 2n 2 Q An 2 uj 2 az 2a 



+ b 2 {a 2 + b 2 y 



\ip\ 2 dx. (A.3) 



It must be that both 2i and I2 are zero if (IA.1I) is to be sat- 
isfied. Supposing that b 7^ then the relationship implied 
by I2 =0, is that 



v 



ul - 2n 2 Q 
a 2 + b 2 



\ip\ 2 dx 



v 



m 2 u 2 az 2a 



(a 2 + b 2 f 



\ip\ 2 dx. 



This result used in the expression for I\ lead to 



h 



\d x ipYdx-k 2 



1 



12 2 



b 2 



\iP\ 2 dx. (A.4) 



v 



However all the quantities in the above expression are 
always positive and, therefore, I\ can never be zero. 
Consequently b must be zero which, in turn, means that 
<7 2 can only be real. 

With 6 = and restoring the definition of a the equation 
I\ — appears now as 







\d x ip\ 2 dx 



2n 2 Q 



4fl 2 UJ 



2, ,2 



(a 2 



\iP\ 2 dx. (A.5) 



Dividing the above expression by \ip\ dx and utilizing 
the definition of Kq given in (|22l) leads to the expressions 
appearing in (|18H20[) . 



